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ABSTRACT 

We compare recent observations of the supernova remnant Gil. 2— 0.3 taken with the VLA during 
2001—02 with images from VLA archives (1984—85) to detect and measure the amount of expansion that 
has occurred during 17 years. The bright, circular outer shell shows a mean expansion of (0.71 ± 0.15)% 
and (0.50±0.17)%, from 20- and 6-cm data, respectively, which corresponds to a rate of 0'.'057±0'.'012/yr 
at 20 cm and 0'.'040 ± 0'.'013/yr at 6 cm. From this result, we estimate the age of the remnant to be 
roughly between 960 and 3400 years old, according to theoretical models of supernova evolution. This 
is highly inconsistent with the 24000 yr characteristic age of PSR J1811— 1925, located at the remnant's 
center, but, rather, is consistent with the time since the historical supernova observed in 386 AD. We 
also predict that Gil. 2— 0.3 is currently in a pre-Sedov evolutionary state, and set constraints on the 
distance to the remnant based on Chandra X-ray spectral results. 

Subject headings: pulsars: individual (PSR J1811 — 1925) — supernovae: individual (Gil. 2— 0.3) — 
supernova remnants 



1. INTRODUCTION 

The process of supernova remnant expansion evolution 
■ has long been the subject of thorough investigation, and as 
a result we can be confident of a few well established facts. 
In the initial free expansion stage, a tremendous amount 
of energy in the form of material ejecta is thrown outwards 
in a (assumed) spherically symmetric explosion, driving a 
shock wave into the ambient medium. Later, when the 
mass of the swept-up material considerably exceeds the 
ejected material mass, the supernova remnant enters the 
Sedov stage (Sedov 1993), as a reverse shock reaches the 
center of the remnant and the forward shock undergoes 
significant deceleration. This simple picture is compli- 
cated by the presence of a pulsar wind nebula expand- 
ing within, but separate from, the SNR bubble (Chevalier 
1984; Reynolds & Chevalier 1984). It eventually encoun- 
ters the reverse shock, which induces complicated radial 
reverberations, before relaxing into Sedov phase expan- 
sion (van der Swaluw et al. 2001; Blondin, Chevalier, & 
Frierson 2001). Each stage of evolution can be described 
by the expansion parameter m, defined by R oc t m , where 
R is the linear radius and t is the remnant age. 

In Gil. 2— 0.3, we find a textbook example of a compos- 
ite supernova remnant (SNR) comprising an extremely cir- 
cular radio and X-ray shell, a pulsar wind nebula (PWN) 
contained within the SNR (Vasisht et al. 1996; Kothes & 
Reich 2001), and an X-ray pulsar PSR J1811-1925 (Torii 
et al. 1997) located at the center. Its high surface bright- 
ness and the extremely centralized position of the pulsar 
within the remnant imply that it is much younger than the 
pulsar's characteristic age, r = 24000 yr (Torii et al. 1999; 
Kaspi et al. 2001), and support the possible association 
with the historical supernova (SN) event of 386 AD (Clark 
& Stephenson 1977). Furthermore, this discrepancy sug- 
gests the pulsar's current spin period is very near its initial 
value and that its spin-down energy, E = 6.4 x 10 36 ergs/s, 
has remained nearly constant since the supernova explo- 
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sion. 

Gil. 2— 0.3 is an ideal SNR for detailed study due to 
the observability of emission from all of its components, 
either at hard X-ray, thermal X-ray or radio frequencies; 
for a description of its radio and X-ray properties see Tam, 
Roberts, & Kaspi (2002) and Roberts et al. (2003). The 
purpose of performing an expansion measurement is to set 
unambiguous upper and lower limits on its age by examin- 
ing the expected behaviour of the outer shock during the 
free expansion (high velocity) and Sedov (low velocity) 
phases. 



2. OBSERVATIONS 

We obtained 20- and 6-cm data of Gil. 2— 0.3 taken dur- 
ing 1984—85 (epoch 1) from the VLA archival database. 
Details of these observations can be found in Table 1 of 
Tam, Roberts, & Kaspi (2002); data from 1985 February 
were omitted from this analysis due to poor calibration. 
Our recent observations made during 2001—02 (epoch 2) 
at 20 and 6 cm (1465 and 4860 MHz, respectively) are 
outlined in Table 1. 

The data processing was performed using standard pro- 
cedures within the MIRIAD package (Sault & Killccn 
1999), in mosaic and multi- frequency synthesis mode. We 
performed calibration and editing on each data set indi- 
vidually, before combining all the data of a particular fre- 
quency band and epoch. The primary gains were deter- 
mined using 3C 286 and 3C 48, and phase calibrations 
were made from observations of 1743—038, 1751—096, 
1751-253 and 1820-254 (J2000). Imaging was performed 
with Robust weighting (Sault & Killeen 1999) as a com- 
promise between maximized signal-to-noise and resolution. 
We utilized the maximum entropy method (MEM) algo- 
rithm for deconvolution (Cornwell, Braun, & Briggs 1999), 
and applied self-calibration iteratively to improve phase 
and amplitude calibrations. 
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3. ANALYSIS 

The VLA and other radio interferometers have been 
used to measure the expansion of many SNR such as 
Gil. 2— 0.3. For a summary of remnants and the tech- 
niques used to study them, see Reynoso et al. (1997) and 
the references therein. Because it is not always possible; to 
make new observations with the same u — v coverage as the 
archival data, it is important to match up as many of the 
properties which might affect the quality of our results 
as best possible, before attempting to directly compare 
the final images from each epoch. We used the final self- 
calibrated epoch 2 clean map, after correcting it for pri- 
mary beam attenuation, as the model for self-calibrating 
epoch 1 data, as described by Masson (1986) in his cross- 
calibration method. The purpose of this step was to apply 
the same residual calibration errors to epoch 1 as existed in 
epoch 2, thereby minimizing the effects of self-calibration 
errors on the final subtracted map. In order to match 
the spatial scales of our images at both epochs, we used 
MIRIAD modelling procedures to spatially filter the u — v 
coverage of epoch 2 data to match that at epoch 1, thus 
creating an epoch 2 dataset with degraded visibility cov- 
erage (Gaensler et al. 1999). 

Rather than subtract our maps in the u — v plane as 
done by Masson (1986), we instead mimicked the direct ap- 
proaches of Moffctt, Goss, & Reynolds (1993) and Reynoso 
et al. (1997) (used to measure the remnant expansions of 
SN 1006 AD and Tycho's SN (3C 10), respectively) who 
adopted procedures outlined by Strom, Goss, & Shaver 
(1982), Tan & Gull (1985) and Dickel et al. (1988). Us- 
ing the MIRIAD task IMDIFF, we fit the five parameters, 
described by the maximum likelihood algorithm of Tan & 
Gull (1985), between our final images: expansion, ampli- 
tude (the mean brightness ratio between epochs), x-shift 
(in the negative right ascension direction) , y-shift (declina- 
tion), and DC offset (the difference in background bright- 
ness levels) . The best-fit geometrical center of the shell was 
found to be < 0'.'5 from the position of PSR J1811-1925. 
We then fixed the amplitude, x-shift, y-shift and offset at 
the fitted values, performed a series of image subtractions 
of epoch 1 from epoch 2, each time artificially scaling the 
epoch 1 image by an expansion factor between 0% and 
1.5% in steps of 0.1%, and examined the radial profiles 
of each difference image for the best-fit expansion factor, 
or more specifically, the profile that most resembled a line 
of zero slope. It was evident when examining the differ- 
ence images that 1.5% was a sufficiently large upper bound 
on the expansion factor. The difference images were con- 
volved with Gaussians whose FWHM were the same as the 
synthesized beams' (19'.'3 x 15'.' 3 at 20 cm, 8'.'3 x 7'/7 at 6 
cm). Figure 1 shows a map of epoch 2 minus epoch 1 at 
20 cm with zero expansion. 

We divided the remnant into 24 wedge-shaped regions of 
15° azimuthally, and found the average flux of each differ- 
ence image in annular ring sections as a function of radius 
between 1' and 3', recalling that the radius of Gil. 2— 0.3 
is ~ 2f25 (Tarn, Roberts, & Kaspi 2002). The purpose 
of averaging over wedges, as opposed to simply taking ra- 
dial cuts, was to smooth out small fluctuations that might 
have corrupted our data (Strom, Goss, & Shaver 1982). 
The next task was to determine which of these profiles of 
average residual flux vs. radius most resembled a flat line. 



A straightforward y 2 analysis was not possible, due to our 
lacking measured uncertainties a on the residual profile 
data points. Therefore, we chose to take a modified ap- 
proach to y 2 -fitting in the hopes of obtaining somewhat 
legitimate error bars. We weighted each profile data point 
according to the total intensity in that region at that ra- 
dius, in lieu of a, and calculated y 2 at each expansion 
value using a flat line as the expected difference profile. 
To find the best expansion and uncertainties for a partic- 
ular region, we divided the \ 2 values by the minimum y 2 
of that region; therefore, when we fitted the values to a 
parabola, its minimum was forced to be near 1. The best 
expansion was determined by the location of the minima, 
and the la errors from the location of y 2 — 2. It should be 
noted that although these error bars do not represent la 
in the traditional sense, they do provide a rough indication 
of each measurement's precision. 

4. RESULTS 

Figure 2 contains a plot of the best expansion rate es- 
timates as a function of azimuthal angle. As expected, 
the largest uncertainies correspond to the regions of the 
shell that are most diffuse, specifically the south-western 
quadrant (between 90° and 180° W of N). We calculate 
the weighted mean expansion rate for the entire shell, and 
find the overall expansion during the ~ 17 yr period sep- 
arating the epochs to be (0.71 ± 0.15)% from 20-cm data 
and (0.50 ±0.17)% from 6-cm data. This corresponds to a 
rate of Of! 057 ± 0'.'012/yr and 07040 ± 0'.'013/yr (from 20- 
and 6-cm data, respectively). The uncertainties represent 
the RMS deviation about the weighted mean. The expan- 
sion parameter m ~ (AR/R)/(At/t) will later be used to 
constrain a rough estimate of the age of Gil. 2— 0.3; here, 
we assume that the age is equivalent to the time since 
SN 386 AD, t = 1616 yr, and calculate m = 0.68 ± 0.14 
(20 cm) and 0.48 ± 0.16 (6 cm). We do not know why 
there exists a general trend in the difference between the 
two sets of results, making the expansion at 20 cm ap- 
pear consistently greater than that at 6 cm; however, we 
note that the 6-cm data is more susceptible to errors in 
primary beam correction and incomplete u — v coverage. 
Even so, the error bars on our data points overlap sig- 
nificantly and the calculated mean expansion values agree 
within our scatter-based uncertainties. 

As an independent check, we compare our measured ex- 
pansion with the general overall expansion for the entire 
remnant as fit by IMDIFF. Our results agree with the fit 
values of 0.7% and 0.5% at 20- and 6-cm, respectively. 

5. DISCUSSION 

In the first phase of SN evolution, the outer shock wave 
freely expands into the surrounding medium such that 
R oc t, and the mass of the ejected material M e j is much 
greater than the mass of swept up material M sw . As 
the forward shock begins to decelerate, a reverse shock 
is driven towards the center of the SNR where it interacts 
with the PWN, if present (Reynolds & Chevalier 1984; 
Blondin, Chevalier, & Fricrson 2001; van der Swaluw et al. 
2001). Eventually these interactions dissipate as the SNR 
enters the Sedov phase. At this stage, M e j « M sw , and 
the expansion is described by the Sedov solution R oc t 0A 
(Sedov 1993). Interactions between the SN shock front and 
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its ambient medium, which has a density profile thought 
to be affected by a circumstellar wind produced by pre- 
supernova mass loss of Type II SN, are currently the sub- 
ject of extensive investigations (Chevalier 1982; Chevalier 
& Fransson 1994; Truelove & McKee 1999). The m values 
predicted for such an environment tend to lie between the 
extreme values of free expansion and the classical Sedov 
solution for a constant density medium. Therefore, we will 
consider these classical definitions for now. 

Assuming that Gil. 2— 0.3 has been in either the free 
expansion or Sedov state continuously since birth, we use 
our above results and the definition of m to calculate upper 
and lower constraints on its age, where AR/R = Ad/ 9 = 
(0.0071 ±0.0015) and (0.0050 ± 0.0017). From our 20-cm 
result, we find upper and lower age limits, corresponding 
to m = 1 and m = 0.4 respectively: ij = 2400 ± 500 yr, 
ti = 960 ± 200 yr. Likewise, from the 6-cm result we find 
i T = 3400 ± 1100 yr and t l = 1400 ± 500 yr. It is immedi- 
ately obvious that the result of this calculation is highly in- 
consistent with the characteristic age of PSR J1811 — 1925, 
t = 24000 yr; it is, however, of the same order of mag- 
nitude as the time passed since SN 386 AD, and indeed 
1616 yr falls easily within these constraints. To eliminate 
any doubt as to the association between the pulsar and 
SNR, see the detailed discussion by Kaspi et al. (2001); 
also discussed are the implications of this age discrepancy 
on young pulsar astronomy. 

Although the PWN is not bright enough to measure 
a rate of expansion in a similar fashion, we find the pres- 
ence of flux outlining the general region of the PWN in the 
zero expansion difference image noteworthy (see Figure 1). 
The positive emission is indicative of positive expansion; 
therefore, we conclude that the reverse shock is not com- 
pressing the PWN shock front and the SNR is not yet in 
the Sedov phase. We look to the study of nonradiative 
SNR evolution by Truelove & McKee (1999), who predict 
analytically and numerically the expected forward blast- 
wave and reverse shock positions throughout a remnant's 
lifetime, for a corrobrative explanation. Based on the fact 
that Gil. 2— 0.3 is likely 1616 yr old, and the expansion pa- 
rameters we find corresponding to that age, we estimate 
roughly that the SNR reverse shock radius is currently 
between 0.5 and 0.8 times the forward blastwave shock ra- 
dius, assuming a typical ejecta mass for Type II SN of 3—5 
M Q . The ratio of PWN to SNR diameter is - 0.28, so the 
reverse shock would not be expected to have reached the 
PWN yet. This agrees with what we observe in Figure 1, 
as well as conclusions outlined in Tarn, Roberts, & Kaspi 
(2002) based on the hydrodynamical simulations of van der 
Swaluw et al. (2001). Furthermore, we refer to Chevalier 
(1982), who predicts that the initial expansion phase of 
a SNR with a red supergiant progenitor can be described 
by a self-similar solution with a value of m — 0.9, as long 
as both the circumstellar material and the stellar envelope 
density distributions are power-laws in radius. This is con- 
siderably higher than our m estimates, based on t ~ 1616 
yr, of roughly 0.48 to 0.68, which suggests that the transi- 
tion from the initial phase to the Sedov phase is well under 
way. 

5.1. Distance Estimate 

The distance to Gil. 2— 0.3 is estimated by consider- 
ing its expansion and angular size. Given the relations 



v = mR/t and R — 8d, where v is the shell's spatial ve- 
locity, is the angular radius, and d is the distance to the 
remnant, it can be seen that 

^"{j-KrTr)- 

To find the velocity of the shell, we consider the Mach num- 
ber of the SNR shock front M — v/c s , where c s = \pypjp 
is the sound velocity (Longair 1994). Here, 7 = 5/3, and 
the ratio of pressure to particle mass density is given by 
p/p = kTi/firrip, where T\ is the temperature of the sur- 
rounding material and fim p is the mean mass per parti- 
cle (/U = 0.6 for cosmic abundances, m p is proton mass) 
(Reynolds et al. 1994). Longair (1997) quotes T 2 /T 1 = 
5M 2 /16 for a strong shock in an ideal gas. We measure 
the temperature behind the shock T\ — 7 x 10 6 K from 
X-ray spectral fits (Roberts et al. 2003); however, T2 is 
the ion temperature and if the electrons are not in full 
thermal equilibrium with the ions, the observed spectrum 
may underestimate the shock temperature, and our dis- 
tance estimate will be too small by a factor of \lTijTx 
(Borkowski, Lyerly, & Reynolds 2001). Combining the 
above information we find 
T 2 5 



Ti 16 \5kTi/3fj,m p y 
which gives a lower bound on the distance estimate: 

_ / mm \ 1/2 f At \ [T£ ( 0.0071 ^ 
\3fimp) {e.AR/Rj~ 6 ]lT x {AR/R) kpC - 

Green et al. (1988) previously estimated a minimum dis- 
tance of ~ 5 kpc to the remnant based on its H I spectrum. 
Fits to the X-ray spectrum with the NPSHOCK model of 
Borkowski, Lyerly, & Reynolds (2001) suggest that the 
electrons are near equilibrium; therefore, the distance de- 
rived assuming total equilibrium should be very close to 
that at near equilibrium, and, hence, not much greater 
than the minimum H I distance. 

6. CONCLUSIONS 

Based on radio interferometric images of SNR 
Gil. 2— 0.3 we have made a simple measurement of the 
outer shell expansion and found a mean rate of 0'.'057 ± 
CC.'012/yr from 20-cm data, and 0'.'040±0 / .'013/yr from 6-cm 
data. If we compare the expected age of Gil. 2— 0.3, deter- 
mined by our measurements, with the characteristic age of 
its associated pulsar PSR J1811 — 1925, we find an order 
of magnitude discrepancy; our result further strengthens 
the growing body of evidence linking Gil. 2— 0.3 with the 
historical SN of 386 AD. The evolutionary status of this 
SNR appears to be pre-Sedov, a conclusion that agrees 
with other observational evidence, as well as theoretical 
arguments. We also estimate the distance to the remnant 
based on its X-ray shock velocity to be > 3 kpc and find 
it consistent with previously published results. 

We wish to thank S. M. Ransom, F. P. Gavriil, V. M. 
Kaspi, M. Lyutikov and S. P. Reynolds for their helpful 
comments and suggestions. The National Radio Astron- 
omy Observatory is a facility of the National Science Foun- 
dation operated under cooperative agreement by Associ- 
ated Universities, Inc. 
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Fig. 1. — 20-cm difference map of epoch 1 image subtracted from epoch 2 image, with no expansion applied. Positive (white) emission 
on the outer shell and negative (black) emission on the inner shell indicate that a noticeable amount of expansion has occured between 
epochs. Contours of the Chandra 0.6—1.65 kcV X-ray image, smoothed with a 5" Gaussian, are shown at levels of 1.5, 2, 3.5, and 5 X 10 — 6 
photons/cm 2 /s/pixcl. The extended structure in the interior of the X-ray shell is thought to be the PWN forward shock; this is discussed in 
Roberts et al. (2003). 
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Fig. 2. — Azimuthal variation of average expansion rate, from 20- (top) and 6-cm (bottom) data. The fitted values (circles) determined by 
a pseudo \ 2 analysis with ~ 1<t error bars are shown (solid lines), as is the weighted mean (dashed line) and its RMS scatter (dotted lines) 
for comparison. 



Table 1 

vla observing parameters for epoch 2 



Observing 


Array 


Frequencies 


Bandwidth 


Time on 


Date 


Config. 


(MHz) 


(MHz) 


Source (min) 


2001 Jun 26 


CnB 


1465 


25 
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2001 Jun 26 


CnB 
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25 
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2001 Aug 03 


C 


1465 


25 
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2001 Aug 03 
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4835, 4885 


25 
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2001 Sep 24 


DnC 
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25 
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2001 Sep 24 
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2002 May 24 


BnA 
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2002 May 24 
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